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For the first time, a next-to-leading BFKL study of the cross section and azimuthal decorrellation 
of Mueller Navelet jets is performed, i.e. including next-to-leading corrections to the Green's 
function as well as next-to-leading corrections to the Mueller Navelet vertices. The obtained re- 
sults for standard observables proposed for studies of Mueller Navelet jets show that both sources 
of corrections are of equal and big importance for final magnitude and final behavior of observ- 
ables, in particular for the LHC kinematics investigated here in detail. The astonishing conclusion 
of our analysis is that the observables obtained within the complete next-lo-leading order BFKL 
framework of the present paper are quite similar to the same observables obtained within next- 
to-leading logarithm DGLAP type treatment. This fact sheds doubts on general belief that the 
studies of Mueller Navelet jets at the LHC will lead to clear discrimination between the BFKL 
and the DGLAP dynamics. 
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1. Introduction 

Soon after the elaboration of QCD, investigations of its behaviour in the high energy regime 
arose. In the semi-hard regime of a scattering process in which s 3> — t, logarithms of the type 
[a v ln(s/|f|)]' 1 have to be resummed, giving the leading logarithmic (LL) Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) [1] Pomeron contribution to the gluon Green's function describing the exchange 
in the ^-channel. The question of testing such effects experimentally then appeared. Based on 
new experimental facilities, characterized by increasing center-of-mass energies and luminosities, 
various observables have been proposed and often tested in inclusive [2], semi-inclusive [3] and ex- 
clusive processes [4]. The basic idea is to select specific observables in order to minimize standard 
collinear logarithmic effects a la DGLAP [5] with respect to the BFKL one. This aims to choose 
them in such a way that the involved transverse scales would be of similar order of magnitude. 

We will here consider one of the most famous testing ground for BFKL physics: the Mueller 
Navelet jets [6] in hadron-hadron colliders, defined as being separated by a large relative rapidity, 
while having two similar transverse energies. One thus expects an almost back-to-back emission in 
a DGLAP scenario, while the allowed emission of partons between these two jets in a BFKL treat- 
ment leads in principle to a larger cross-section, without azimuthal correlation between them. The 
predicted power like rise of the cross section with increasing energy has been observed at the Teva- 
tron /^-collider [7], but the measurements revealed an even stronger rise than predicted by BFKL 
calculations. Besides, the leading logarithmic approximation [8] overestimates this decorrelation 
by far. Improvements have been obtained by taking into account some corrections of higher order 
like the running of the coupling [9]. Calculation with the full NLL BFKL Green's function [10] 
have been published recently [11]. We here review results of Ref. [12] on the full NLL BFKL 
calculation where also the NLL result for the Mueller Navelet vertices [13] is taken into account. 



2. NLL calculation 

The kinematic setup is schematically shown in Fig. [I]. The two hadrons collide at a center 
of mass energy s producing two very forward jets, whose transverse momenta are labeled by Eu- 
clidean two dimensional vectors k/j and k/2, while their azimuthal angles are noted as <j)j t \ and 
07,2- We will denote the rapidities of the jets by yj \ and yj i2 which are related to the longitudinal 
momentum fractions of the jets via xj = ^e yj ■ Since at the LHC the binning in rapidity and in 
transverse momentum will be quite narrow [14], we consider the case of fixed rapidities and trans- 
verse momenta. Due to the large longitudinal momentum fractions xj\ and xj 2 of the forward jets, 
collinear factorization holds and the differential cross section can be written as 

.,- I./ . . =1 fd*! f l dx 2 M Xl )Mx 2 ) dd » b , (2.1) 

d | ky ; 1 1 d | k^ 2 1 dyy, i dy J;2 ~b Jo Jq d | ky, i | d | k^ )2 1 dyy, i d.y^ j2 

where f a t, are the parton distribution functions (PDFs) of a parton a (b) in the according proton 
(which are renormalization scale ju« and factorization scale jip dependent). The resummation of 
logarithmically enhanced contributions are included through kj -factorization: 

,,, |,| 1 dg ?, ~ A = / #/,id<^, 2 [ d 2 k 1 d 2 k 2 V a (-k 1 ^ 1 )G(k 1 ,k 2 ,f)y b (k 2 ^ 2 ), (2.2) 

where the BFKL Green's function G depends on s = x\X 2 s. The jet vertices V a ^ were calculated at 
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Figure 1: Schematical illustration of the kinematics. 



NLL order in Ref. [13]. Combining the PDFs with the jet vertices one writes 

da 

d|k /a |d|k 7i 2|d^idj7,2 

where $(k, j2) x 7)2) k 2 ) = /dx 2/ (x 2 )V(k 2 ,x 2 ). (2-3) 



J d0 7il d0 7 , 2 J d 2 k 1 d 2 k 2 0(k / , 1 ,x 7 , 1 ,-k 1 )G(k 1 ,k 2 ,f)<I>(k 7i2 ,x / , 2 ,k 2 ) 



In view of the azimuthal decorrelation we want to investigate, we define the following coefficients: 



=Jd(()j A d0 7 , 2 cos (m((j)j A - (j) Ja - n)) 7d 2 ki d 2 k 2 <J>(k 7! i,xy ; i, -ki)G(k 1 ,k 2 ,f)<I>(k /j2 ,X7 j2 ,k 2 ), 
from which one can easily obtain the differential cross section and azimuthal decorrelation as 

— — d ° = % and (cos(m<p)) = (cos (m(<j>ji - (j> J)2 ~^))) = (2.4) 

The guiding principle is then to rely on the LL-BFKL eigenfunctions 

^, v (k 1 ) = -L(k 2 )' V ^^', (2.5) 

7l\/2 

although they strictly speaking do not diagonalize the NLL BFKL kernel. In the LL approximation, 

(o(my) 



where 



r ( § \ W(m ' v > 

tf m = (4-38 mfi ) ydvC^vG^il^ijC^vO^l^f-J , (2.6) 

C m>v (\kj\,xj) = J d0 / d 2 kdx/(x)y(k,x)£ m , v (k)cos(m0 7 ), (2.7) 

and a>(n,v) = a s Xo(\n\,j + iv), with Xo(n,f) = 2¥(1) - ¥ (7+ §) - ¥ (l - 7+ §) and a,. = 
N C CC S / 71. The master formulae of the LL calculation ( |2.6| , |2.7| ) will also be used for the NLL cal- 
culation. The price to pay in order that E n V remains eigenfunction is to accept that the eigenvalue 
become an operator containing a derivative with respect to v. In combination with the impact fac- 
tors the derivative acts on the impact factors and effectively leads to a contribution to the eigenvalue 
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which depends on the impact factors. 

At NLL, the jet vertices are intimately dependent on the jet algorithm [13]. We here use the 
cone algorithm, which is expected to the used by the CMS collaboration. At NLL, one should also 
pay attention to the choice of scale sq. We find the choice of scale sq = \Ao,i *o,2 with so l = tr^kj 1 
rather natural, since it does not depend on the momenta ki 2 to be integrated out. Besides, the 
dependence with respect to so of the whole amplitude can be studied, when taking account the fact 
that both the NLL BFKL Green function and the vertex functions are so dependent. In order to study 
the effect of possible collinear improvement [15, 16], we have, in a separate study, implemented 
for n = the scheme 3 of Ref. [15]. This is only required by the Green function since we could 
show by a numerical study that the jet vertices are free of 7 poles and thus do not call for any 



collinear improvement. In practice, the use of Eqs. ( |2.6| , |2.7[ ) leads to the possibility to calculate 
for a limited number of m the coefficients C m , v as universal grids in v, instead of using a two- 
dimensional grid in k space. We use MSTW 2008 PDFs [17] and a two-loop strong coupling 
with a scale pL R = \/|k/~i | • |k/.2| . Although a BFKL treatment does not require any asymmetry 
between the two emmited jets, in order to compare with DGLAP NLO approaches [18] obtained 
through the NLL-DGLAP partonic generator Duet [19], for which symmetric configurations lead 
to instabilities, we here display our results for |k/j| = 35GeV, |kj^| = 50GeV, and compare them 
with the DGLAP NLO prediction (see Ref. [12] for symetric configurations). 



3. Results 

We now present our results for the LHC at the design center of mass energy y/s = 14TeV. 
Motivated by a recent CMS study [14] we restrict the rapidities of the Mueller Navelet jets to the 
region 3 < \yj\ < 5, thus limiting the relative rapidity Y between 6 and 10 units. Fig. ||a and ||b 
respectively display the cross-section and the azimuthal correlation 1 . This explicitely shows the 
dramatic effect of the NLL vertex corrections, of the same order as the one for the Green function. 
In particular, the decorrelation based on our full NLL analysis is very small, similar to the one 
based on NLO DGLAP. The main source of uncertainties is due to the renormalization scale }JLr 




Figure 2: Differential cross section and azimuthal correlation in dependence on Y for |k/ i| = 35GeV, 
|k/ 2| = 50GeV. The errors due to the Monte Carlo integration are given as error bands. As dots are shown 
the results of Ref. [18] obtained with DlJET [19]. 

and to the energy scale y/sfj. This is particularly important for the azimuthal correlation, which, 

'We use the same color coding for both plots, namely blue shows the pure LL result, brown the pure NLL result, 
green the combination of LL vertices with the collinear improved NLL Green's function, red the full NLO vertices with 
the collinear improved NLL Green's function. 
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when including a collinear improved Green's function, may exceed 1 for small [Ir = jJ.p. 

In conclusion, contrarily to the general belief, the study of Mueller Navelet jets is probably 
not the best place to exhibit differences between BFKL and DGLAP dynamics. 

Work supported in part by the Polish Grant N202 249235, the grant ANR-06-JCJC-0084 and 
by a PRIN grant (MIUR, Italy). 
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